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1. Introduction 



Non-Gaussianity has become a very important probe into the physics in the early 
universe. A well-understood non-Gaussianity has a local form which says that the curvature 
perturbation can be expanded to the non-linear orders at the same spatial point 



9 



C(x) = C s (x) + -/jvlC(x) + — ffiVLC(x) + • 



where f^L and gNL are the non-Gaussianity parameters which set the sizes of bispectrum 
and trispectrum respectively. Single-field inflation predicts /nl ~ 0(n s — 1) which is con- 
strained to be much less than unity. A convincing detection of local form non-Gaussianity 
will rule out all single-field inflation models. 

A large local form non-Gaussianity can be generated by the isocurvature field(s) at 
the end of multi-field inflation 0, ||, Q or deep in the radiation-dominant era, such as 
curvaton model H ^, ^, ||, [l0|, [ll], [l^, [0| 14] . In the literatures, the non-Gaussianity 
parameters are assumed to be scale-independent. However, recently ones found that a 



scale-independent Jnl is not a generic prediction of inflation (15, p], 17, 18, 19, 20, 21]. 



For simplicity, the scale-dependent Jnl and g^L are parameterized as follows 

/ k \ n fNL + l a fNL ln ^ 

/atlO) = fNL(k p ) [ -r- 



ln- 



gNh(k) = g NL {kp) ( — 



fl.2l 



Cl-3) 



where k p is a pivot scale. If such a scale dependence is not too small, it can be possibly 



detected by the forthcoming experiments. For example, in [23], the authors showed that 



lr The scale dependence of f'ST' atera was discussed in 
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Planck [ P4 | and CMBPol [25] are able to provide a 1-a uncertainty on the spectral index 
of fjsri for local form bispectrum: 

50 1 



An 



0.1 



fNL sf fsky 



and 



An fNL ~ 0.05- 



50 



for Planck, 



for CMBPol, 



(1.4) 



(1.5) 



fNL \J fsky 

where f s k y is the sky fraction. The effects on the large-scale structure from scale-dependent 
fNL were discussed in [26, 27]. From the theoretical point of view, the leading order of 
non-Gaussianity can be higher terms, such as that set by gNL-, not fNL- Assuming f^L = 0, 
the current observational limit for g^L are —7.4 x 10 5 < g^L < 8.2 x 10 5 (95 % C.L.) from 
cosmic microwave background observations [^] and —3.5 x 10 5 < g^L < 8.2 x 10 5 (95 % 
C.L.) from large scale structure observations [29]. Planck will reduce the uncertainty of 
gNL to A<7tvl = 1-3 x 10 5 [29]. More studies on fingerprints of the scale-dependent g^L in 
CMB and large-scale structure are called for in the future. 

A large non-Gaussianity implies a strong interaction between different cuvature per- 
turbation modes. However, it can be generated by the isocurvature field without self- 
interaction at all. In order to dig out more physics about the isocurvature field, we need 
some new observables. In [16, 17, 18, 19 1, we found that f^L generated by a free isocur- 



vature field is scale independent and the spectral index and running of /jvl are good 
discriminators to the self-interaction of isocurvature field. We extend our discussions on 
the scale dependence of /jvx in [19 to g^L in this paper. 

Our paper is organized as follows. In Sec. 2, we use the method in [19| to derive the 
spectral and running of gNL from a single isocurvature field. In Sec. 3, for an example, we 
apply our formula to investigate the curvaton model with near quadratic potential. More 
discussions are contained in Sec. 4. 



2. The spectral index and running of g ml from an isocurvature scalar field 

In this paper we consider that the curvature perturbation is generated by the quantum 
fluctuation of an isocurvature field a which slowly rolls down its potential during inflation. 
Its dynamics is governed by 

(2.1) 

3# V ' 

where Via) = dV(a)/da and H is the Hubble parameter. 

The gravitational dynamics itself introduces important non-linearities, which will con- 
tribute to the final non-Gaussianity in the large-scale CMB anisotropics. Based on the 
so-called 5N formalism pOfl , the curvature perturbation produced by the isocurvature field 
can be expanded to the non-linear orders as follows 

C(*/,x) = A^i/^M^xH^A^/,^)^ , (2.2) 
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where N >a , N >aa and N^ aa are the first, second and third order derivatives of the number 
of e- folds with respect to a respectively. Here tf denotes a final uniform energy density 
hypersurface and ij labels any spatially flat hypersurface after the horizon exit of a given 
mode. Similar to [19], tj is set to be t*(k) which is determined by k = a(t*)H* for a given 
mode with comoving wavenumber k. Therefore the amplitude of the curvature perturbation 
is given by 



A 2 



HA' 
2vr J 



(2.3) 



Gravitational waves were also generated during inflation and the amplitude of its power 
spectrum takes the form 

Hi 



A 



T 



7T 2 /2' 



(2.4) 



Here we work on the unit of M p = 1. The scale dependence of gravitational wave pertur- 
bation is measured by i%t which is defined by 



tit 



where 



dA 2 T 
dink 



-2e H , 



(2.5) 



(2.6) 



Usually we introduce a new quantity, the tensor-scalar ratio r^, to measure the amplitude 
of gravitational waves: 



r T = A^/A 



8 



n 



From Eq. 



and 



1 2j), the non-Gaussianity parameters are given by 

- 6 N 2^q ■ 



9NL 



25N aaa (U 



(2.7) 



(2- 



(2.9) 



54 N%(Q • 

For working out the scale dependence of gNL, we follow the method in [19] and intro- 
duce a new time i r (> i*) which is chosen as a time soon after all the modes of interest exit 
the horizon during inflation. The value of a at t r is related to that at time by 

rtr r H 

(2.10) 



da 



V a) 



3H{t) ' 



Therefore we have 



da r 




V'(a r ) 


da* 


u 


V>(a*) 


da r 




V'(a r ) 


dt* 




3H(U) 



(2.11) 
(2.12) 
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Considering 



d F ( Jr ) _ dF(a r ) da r | da r 
dt* da r a* dt* 



and = —V'(a^), one finds 



d :n°r) = -eAtF^) = 0, (2.14) 



din A; H*dt* 

which implies that F(a r ) is scale independent. Taking into account that a r is a function 
of cx*, we have 

_ &a r dN(cr r ) o da r d 2 a r d 2 N(a r ) ( da r \ 3 d 3 N(a r ) 
da 3 8a r + daedal da? + \daj da? ' 

Since dN Q° r ^ , 9 gffi 1 ^ an d 8 gffi^ are scale independent, one obtains 
dlnN :t7 (U 



dink 



Vaa, (2-18) 



f»5^_^ + *fc£> (2 . 19) 

d\nN aar7 (Q _ N aa (U) N a (Q 

— — - 6i] aa + 3773— — + £ 4 — — , (2.20) 

dink N, aaa (U) N jaaa (U) 

where the slow-roll equation for a is adopted and 

^ = ^Hf' % = ^T' e4_ ^r- (2 - 21) 

From the above results, the spectral index of A^, /nl and g^z, are respectively given by 

din At 

n s = l + = 1 + 2rj aa - 2e H , (2.22) 



_ dln\f NL \ N :a (Q 

= -^nTk- = (2 - 23) 

and 

_ dln\g NL \ N aa (Q N a (Q 

n ^ = ~d^r = 3 ^ + (2 - 24) 

= 2 V Ln ^ + W2^- (2 - 25) 

9NL 462 g NL 
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Our results are the same as those in [p^7|| . The scale dependence of the non-Gaussianity 
parameters from an isocurvature field is proportional to third and fourth order deriva- 
tives of its potential. Therefore the spectral indices of Jnl and g^L are really the good 
discriminators to measure the self-interaction of such an isocurvature scalar field. 

In fact the indicies n s and ut may be scale dependent as well. Their scale dependences 
are measured by the so-called runnings which are defined by 



dn s 
dink' 



and ax 



driT 
din A; 



(2.26) 



The dynamics of slow-roll inflation is governed by the inflaton field (j) whose potential is 
denoted as V/(</>) and then we have 



where 



Considering 



we obtain 



den 
din k 



-2€HVU + 4e tf' 



drj a 



dink 



Vi(4>) 



2e H r]aa -£3, 



(2.27) 



(2.28) 



(2.29) 



a s = a T + 4:e H r]aa - 2£ 3 . 



(2.30) 
(2.31) 



where 



V'V" 



(2.32) 



In the future, the spectral index and running of curvature perturbation might be measured 
precisely. There might be an opportunity to measure by CMBPol |25[ as long as is 



not too small [32|. One can imagine that it is quite hard to measure ax- Prom the above 
formula, rj^ makes contribution to a s as well. So it is difficult to re-construct the potential 
of isocurvature field from n s and a s . 

In [jl!]] the running of spectral index nj NL from an isocurvature field is derived as 
follows 



dink 



(2eu - rjaa - 77 4 )n/ A 



where 



V4 



3H 2 V" 



(2.33) 



(2.34) 
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Similarly, we define the running of n gNL , namely 



Taking into account 



and 



the running of g^L becomes 



dn 



a 



9nl 



9nl 



dink 



dink 



drx 
dink 



(2e H ~ 



-2r T r)aa, 



(2.35) 



(2.36) 



(2.37) 



a 



9NL 



n gNL (2r] aa - 2e H ) - n 



2 

9NL 



+ 



9NL 



where 



25 1 t 



(2.38) 



(2.39) 



If n 9NL ~ 0(1), oig NL is roughly the same order of n 9NL , and hence n gNL is not a reliable 
quantity to characterize the scale dependence of qnl an y more. Actually there is one 
another non-Gaussianity parameter tnl which measures a special local form trispectrum. 
In our setup, tjvl is not an independent parameter and it is related to /jvl by tjvl = 



(I/atl) 2 . Therefore 



n 



TNL 



2rif NL and a 



TNL 



2a fNL ■ 



Combining with Eq. (|2.22|), the running of spectral index of g^L becomes 



a 9NL — n 9iVi( n s l + n 3jvz,) 

+ — — 2f NL n fNL (rjtra -m + n fNL ) 



9NL 



25 
432 



(2.40) 



For the special case with /jvl = 0, the spectral index and running of g^L are simplified to 
be 



n 



9NL 



25_ r T 
^2 U g NL 



a 



9nl 



n 



9nl 



(1 - n s - £5 - n 



9nl 1 



(2.41) 
(2.42) 



Now for an isocurvature field with renormalizable potential V(cr) ~ a n , where n is an 
integer not larger than 4, £5 = and we obtain a consistency relation 



a 



9NL 



n 



9nl 



(1 - n s 



n 



9nl j 



(2.43) 



We hope that the running of g^L can be detected as well even though it is expected to be 
very difficult in the near future. 
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For an instance, we consider an isocurvature field a who has a polynomial potential, 

V(a) = \m 2 a 2 + Xm 4 (-Y . (2.44) 



2 \m/ 

Similar to ©0, (jj, we introduce a new parameter to measure its self-interaction term 
compared to its mass term as follows 



n-2 



s = 2\ — . (2.45) 



Here s > — 2/n. Otherwise the isocurvature field will run away to the infinity. Accordingly, 
tyro-) ??4 and £4 are given by 

n(n — 1) 
1 + ^— 

and 



Vera = 11mm ( H S ) (2.46) 



n(n- l)(n-2) s n(n - l)(n - 2)(n - 3) s 

% = ~ Vmm — , £4 = Vmrn^, (2.47) 

2a* 2 ai 

where r\ mm = m 2 /3H' 2 . Here m 2 is assumed to be positive. The spectral indices of /jvl 
and g^L are 

/at s 5n(n — l)(n — 2) H* s ,„ lnl 

nfsi. = Sign JV !(T V 2 -Vmm — (2-48) 

. . 5n(n - l)(n - 2) H*f NL 
n gNL = sign(iV CT ) *nm-— ' 

I 25»(n-l)(n-2)(n-3) 

< H *J^ i the second term on the right hand side of Eq. ( [2.49D becomes 

dominant. Since both rif NL and n gNL are proportional to r] mm , the scale dependences of 
/nl and g^L are detectable only when the mass of isocurvature field is not too small 
compared to the Hubble parameter during inflation. 

Because the Compton wavelength of isocurvature field is large compared to the Hubble 
size during inflation, its quantum fluctuations leads to a typical vacuum expectation value 
of a ||, i.e. 



/ 3 H 
o* = \ L 2n , w r — ■ (2.50) 
y 87r 2 (l + §s) m 

Taking the above result into account, nj NL and n gNL are simplified to be 

n fNL = sign(iV )<T ) • 2.3 x 10 3 n(n - l)(n - 2)^- (^-\ s^l + ^s, (2.51) 



n gNL = sign(7V ff ) • 4.6 x 10 J n(n - l)(n - 2)^ — sJl + -s 

9NL \H*J V 2 

+ 2.1 x 10 7 n(n-l)(n-2)(n-3)— (-^) s(l + (2.52) 

Sjvl \H*J 2 
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Here we adpot WMAP normalization: = 4.96 x 10 5 [31]. For a model with detectable 
n fNL by PLANCK, |n /jVL • f NL \ > 5 |23| which implies 



m 



> 0.13 



n{n — l)(n — 2)|s| 



n 

1+ 2 5 



-1/3 



(2.53) 



For the model with sign(A r j(T ) > 0, the region of n 9NL -gNL with detectable nj NL is illustrated 
in Fig. [l[ Because a is considered to be a general isocurvature field, one cannot figure out 
the value of /jvx only from m/H* and s. The current constraints on f^L from WMAP [31] 
is Jnl = 32 ± 21. So, for example, we consider /jvx = 10 and 50 which correspond to the 
black and blue curves in Fig. [j], respectively. We find that the lower bound on \n 9NL ■ gNh\ 
is not sensitive to the value of Jnl- On the other hand, the mass of isocurvature field 
is assumed to be smaller than the Hubble parameter during inflation and then Eq. ( 2.53j ) 
implies that |s| cannot be too small. 

3. n 9NL in the curvaton model with near quadratic potential 

We start with the curvaton whose potential takes the form in Eq. (|2.44| ). In this section, 
we focus on the case in which the self-interaction term is much smaller than the mass term. 
For \s\ <C 2/n and m 2 > 0, the amplitude of scalar power spectrum, /jyx and gNL are 
respectively given by 



2 H*/2n 
o o * 



/ 



NL 



4/d 
25 



(l + h 2 



5/d 
6 



9NL= 54 



4/B 



(h 3 + 3h 2 ) - -^-(1 + /i 2 ) + J(l - 9/i 2 ) + 10/ D + 3/|, 



where fo = t~^z , and 



9 
/;2 

/;3 



_ w(xp) + n(n - l)g(n, x )s/2 
w(x ) + ng(n,x )s/2 
-w(xp) + ng(n,x )s/2 
(w(x ) + n(n - l)#(n, x )s/2) 2 
(w(x ) + ng(n,x )s/2) 2 



n(n — l)(n — 2)g(n, xq)s/2, 
n(n — l)(n — 2)(n — 3)g(n, xq)s/2. 



(3.1) 
(3.2) 
(3.3) 

(3.4) 
(3.5) 
(3.6) 



(w(cco) + n(n - l)g(n, x )s/2) 3 
See, for example, [[ll|] in detail. Here Q a ,D is the fraction of curvaton energy density in the 
total energy density budget at the time of its decay, 

w(x ) = 2 1 / 4 r(5/4)xo 1/4 J 1/4 (x ), (3.7) 

and 



Jl/4(xo) J J^ 4 \x)Y 1/4 {x)x^ n y i dx 
Yyiixo) r°J!)4^)^ (6 " n)/4 H , 



(3.8) 
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Figure 1: The regions of n 9NL ■ @nl above the curves correspond to the cases with a detectable 
nf vi- 



and xq = mto = 1 denotes the time when curvaton starts to oscillate. Taking Eq. (| 



into account, the tensor-scalar ratio becomes 

For a sub-Planckian value of <r*, the tensor-scalar ratio is much smaller than one if fo is 
not too small. 

In order to obtain a large non-Gaussianity, fu should be smaller than one. For the 
curvaton model with near quadratic potential, the non-Gaussianity parameters are approx- 
imately given by 

f NL ~ jL(l + h 2 ), (3.10) 
25 

9NL k ^772-^3+3/12). (3.H) 

One can easily calculate the spectral indices of Jnl and gNL- 

n(n — l)(n — 2) s 
n fNL * ~ 2 r l m m q{1 + h2 y (3.12) 

q (l + hv) + ~-l\. (3.13) 



71 ~ 



n(n — l)(n — 2) 3s 



q 2 {h 3 +3h 2 ) 



The spectral indices of both Jnl and ^tvl are independent on fjj. Here ray^ and n 9NL 
are illustrated in Fig. |2| for s < and s > respectively. In this setup there are four 
model-dependent parameters: n, s, rj mm and They can be fixed by four observables: 
/iVXj n /jvLi an d n gjvL- The spectral indices of Jnl and gNL are very useful for us to 
re-construct curvaton potential. 

We also notice that /nl can be tuned to be zero even when fu <C 1 for the curvaton 
model with near quadratic potential. In this case the higher order terms, such as gNL, 
still leads to a large non-Gaussianity which can be potentially detected by the forthcoming 
observations as well. Now h 2 = —1, but gNL can still be quite large 

25 /13-3 

9NL - 24^- (3-14) 

The spectral index of g^L is simplified to be 

1 s 
u 9nl = ~2 n ( n ~ J )( n ~ 2 )( n ~ 3)r? mm ^ 2( , 3 _ (3.15) 

Here, for a given n, s can be fixed by the condition h 2 = — 1. Our numerical results are 
summarized in Fig. ||. From Eq. ( 3.15| ), n gNL is proportional to r/ mm . Larger r) mm , larger 



n gNL . On the other hand, if and tit are not too small, they might be measured in the 
near future [32], and then rj acr and rj mm can be fixed by 770-0- = (n s — 1 — tit)/2 and Eq. ( 2.46|) . 



The order of magnitude of both r\ aa and r\ mm is expected to be 0(1O -2 ). Therefore the 
spectral index of gNL in the curvaton model with /nl = is negative and its order of 
magnitude is roughly — 0(1O~ 2 ). 
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Figure 2: The values of nf NL /rj mm and n gNL /r] mm in the curvaton model with a polynomial 
potential. The solid, dashed and dotted lines correspond to n = 4, 6, 8 respectively. 



4. Discussions 

Inflation is driven by the vacuum energy of inflaton field. The vacuum energy density is 
almost a constant with the expansion of universe and then the Hubble parameter almost 
did not decrease during inflation. However it is not an exact constant, but slowly decreased. 
WMAP data implies that the power spectrum of curvature perturbation is just near scale 
invariant (n s = 0.968 ± 0.012), not exactly scale invariant. The tilt of power spectrum 
comes from the fact that the inflaton potential is not exactly flat. 

A large local form non-Gaussianity can be naturally generated by an isocurvature field 
on the super-horizon scales. If the isocurvature field is a free field (without self-interaction), 
the non-Gaussianity parameters, such as /nl and gNL, should be scale independent. More 
generally, one may expect that the isocurvature field self-interacts with itself and then the 
non-Gaussianity parameters also depend on the wavelengths of perturbation modes. On 
the other hand, one can learn how the isocurvature field interacts with itself from the scale 
dependences of the non-Gaussianity parameters. 

In this paper we derive the spectral index and running of g^L from a general isocur- 
vature field. The typical region of n gNL for the model with detectable nf NL is illustrated 
in Fig. [y. Applying our results to the curvaton model with near quadratic potential, we 
find that both nj NL and n gNL are independent on Therefore one can usually tune the 
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Figure 3: The values of s, rj aa /rj : 



mm, i'j N [ 



/i]mm and gNL*fr> in the curvaton model with /jvl = 0. 



free parameter /o to achieve a large non-Gaussianity. However, in curvaton model, one 
can tune Jnl to be zero even when fjj <C 1, but g^L ~ V/d is still large. In this special 
case, n ffJVL ~ — 0(1O~ 2 ). 

In the literatures, the near scale- invariant variables are usually re-parametrized by 



o{k) = o(k p ) [ A 



\aln- 



(4.1) 



This expression is reliable around k 
follows 



k p . Actually one can re- write the above formula as 



0(k) = 0(k p 



. k 1/ k 
1 + n in 1 — a I In — 

kp 2 \^ kp 



+ ... 



(4.2) 



where a = n 2 + a is the new "running" of spectral index n. If we adopt this new definition, 
the terms — n? and —n gNL on the right hand sides of a/ JVL and a SJVL in Eqs. ( 2.33| ) 



9nl 



and ( p.38 ) can be absorbed into a/ JVL and a 9jVL respectively. An advantage of this new 



definition is that a.f NL and a 9NL are respectively much smaller than nf NL and n gNL even 
when rif NL , n gNL ~ 0(1). 

Once the large local form non-Gaussianity is confirmed by the forthcoming observa- 
tions, the scale dependence of the non-Gaussianity parameters will be the next important 
issue. How to search for the signal of scale dependence of g^i m the CMB and large-scale 
structure data is still needed to be done in the near future. 
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